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ABSTRACT 



Context. Evolutionary synthesis models are the prime method to construct spectrophotometric models of stellar popu- 
lations, and to derive physical parameters from observations through comparison with such models. One of the basic 
assumptions for evolutionary synthesis models so far has been the time-independence of the stellar mass function (ex- 
cept for the successive removal of high-mass stars due to stellar evolution). However, dynamical simulations of star 
clusters in tidal fields have shown the mass function to change due to the preferential removal of low-mass stars from 
clusters. 

Aims. Here we combine the results from dynamical simulations of star clusters in tidal fields (and especially the derived 
parametrisation of the changing mass function) with our evolutionary synthesis code GALEV to extend the models by 
a new dimension: the total cluster disruption time. 

Methods. Following up on our earlier work, which was based on simplifying assumptions of the mass function time 
dependence, we reanalyse the mass function evolution found in N-body simulations of star clusters in tidal fields, 
parametrise it as a function of age and total disruption time of the cluster and use this parametrisation to compute 
GALEV models as a function of age, metallicity and the total cluster disruption time. 

Results. We study the impact of cluster dissolution on the colour (generally, they become redder) and magnitude (they 
become fainter) evolution of star clusters, their mass-to-light ratios (which can deviate by a factor of ~2 - 4 from 
predictions of standard models without cluster dissolution), and quantify the effect on the cluster age determination 
from integrated photometry (in most cases, clusters appear to be older than they are. Depending on the filter set 
available and the evolutionary stage of the cluster, the age difference can range from 20 to 200%). By comparing our 
model results with observed M/L ratios for old compact objects in the mass range we find a strong 

discrepancy for objects more massive than lO'' M0, in a sense that observed M/L ratios are higher than predicted by 
our models. This could be either caused by differences in the underlying stellar mass function or be an indication for 
the presence of dark matter in these objects. Less massive objects are well represented by the models. 
The models for a range of total cluster disruption times and metallicities are available online, at 
http://www.phys.uu.nl/~anders/data/SSP_varMF/ and http://data.galev.org, and will be made available via CDS. 
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1. Introduction 



Since the pioneering work by Tin sley (jTinslevI Il968t 
iTinslev fc GunnI 1X9761 : iTinslevI [T980( ). evolutionary syn- 
thesis modelling has become the method-of-choice to 
predict spectrophotometric properties of stellar popula- 
tions. Currently popu lar models include STARBURStQQ 
(iLeitherer et al.l I1999D. GALAXEV (Bruzual fc Chariot 



20031). GALEV (I Anders fc Fritze-v. AlvenslebenI 120031: 



Bicker et al] |2004 . pegase (iFioc fc Rocca-Volmerangd 



19971 ). and the Maraston models ( Marastonl2005[ ). all giving 
predictions for single-age populations , so-called "Simple 
Stellar Populations" (SSPs). In addition, GALAXEV, 
PEGASE and GALEV also provide models for populations 
with arbitrary extended star formation histories (SFH, like 
galaxies), while STARBURSt99 only allows for an extended 
constant SFH. Comparing predictions from these models 
with observations allows to derive basic physical param- 
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eters of the stud ied system (e.g ., among many others , 
[ Bicker et al. 200 21 iKassin et alj[200 3: Ande rs et al.ll2004bt 
Ide Griis et aL,,20qlTK'undu et al.ir2005; ,de Griis fc Anderd 



20061 : ISmith et al.ll2007l) . 



While the specific input physics (e.g. the choice of stellar 
isochrones and spectral libraries, the inclusion of gaseous 
emission) and implementation varies among the models, 
some basic techniques and limitations are inherent to all of 
them: assigning a spectrum to each star along the isochrone, 
weighting them according to a chosen stellar initial mass 
function (IMF) and integrating along the isochrone (and 
over the SFH, if applicable) results in the integrated prop- 
erties of the stellar population at a given age. For all cur- 
rently available models, the stellar mass function (MF) is 
time-independently fixed at its initial value, the IMF. 
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Cluster disruptior(3 has become a well-studied phe- 
nomenon. It can be observed both in the earliest phases 
in a cluster's life (the so-called "infant mortality" caused 
by the removal of gas left over from the cluster forma- 
tion pr ocess by stellar winds and/or the first supern ovae, 
see e.g. iLada fc Ladall2003l: iBastian fc Goodwinll2006l ) and 
for old clusters (e.g. the prominent tidal tails of the 
Milky Way globular cluster Palomar 5. lOdenkirchen et al.l 
l2003ii . Age and mass distributions of a whole star clus- 
ter system can be used to determine the typical dis- 
ruption time of clusters of a given mass in this cluster - 
system (Boutlou kos fc LamersI [20031 iLamers et all l2005bt 
iGieles et al. 2005r This cluster disruption time is predom- 
inantly determin ed by the external tid al field the cluster is 
experiencing (see ILamers et al" 2005 b^). the local density of 
giant molecular clouds (iGieles et al. 2006) and the occur- 
rence of spiral arm passages (jGieles et al.l [2007! ). In addi- 
tion, the cluster loses mass due to stellar evolution. While in 
the case of "infant mortality" the cluster is likely (almost) 
completely disrupted (although a bound core might remain, 
see e.g. Bast ian fc Go odwin 2006 for "infant weight loss"), 
cluster dissolution in a smooth external tidal field is a more 
gradual process and accompanied by perpetual dynamical 
readjustment within the cluster. The latter is characterised 
by a mass-dependent probability to remove a star from a 
cluster: due to energy equipartition massive stars tend to 
sink towards the cluster center, while low-mass stars are 
driven outwards where they are rnore easily removed by 
the surrounding tidal field ([Henonl 1X9691 : ISpitzer fc Shulll 
119751: iGiersz fc Heggielll997[ ). The resulting radial depen- 
dence of the mean stellar mass inside a cluster is called 
"mass segregation". Mass segregation established by the 
very star formation process itself is referred to as "primor- 
dial mass segregation" (for observa tional evidence of "pri- 
mordial mass seg regation" see e.g. iGouliermis et al.ll2004t 

IChen et al.ll2007h. 

iBaumgardt fc Makinol (|2003D (from now on: BM03) per- 
formed the first (and, so far, most extensive) quantitative 
large-scale study of how the stellar MF inside a star cluster 
changes due to dynamical cluster evolution in a tidal field. 
They confirmed earlier findings for a preferential loss of low- 
mass stars and derived a formula describing the change in 
MF slope for low-mass stars. However, their derived formula 
(formula (13) in BM03) only applies to stars with masses ^ 
0.5 M0, while the effect is pronounced also for higher-mass 
stars (see BM03 Fig. 7), which dominate the flux emerg- 
ing from the cluster (for ages shorter than a Hubble time). 
BM03 performed their simulations for clusters not pri- 
mordi ally mass-segregated. Recently, in IBaumgardt et al.l 
(|2008D . they studied also the dissolution of initially mass- 
segregated clusters (with a simplified initial setup differ- 
ent from the BM03 simulations, hence we can not combine 
these sets of simulations), finding an ev en stronger MF evo- 



studied the evolu- 



tnese sets oi simula tions), nnamg an eve 
lution than BM03. iMarks et al.l (|2008D 
tion of the stellar MF inside star clusters during the gas 
removal/ "infant weight loss" phase, and found it to also 
preferentially remove low-mass stars, leading to a flatten- 



^ With disruption we encompass all kinds of different cluster 
mass loss and destruction events (e.g. single disruptive encoun- 
ters with giant molecular clouds, infant mortality, final cluster 
"death"). Dissolution stands for any gradual destruction pro- 
cess, e.g. mass lost due to stellar evolution, tidal dissolution or 
multiple weak encounters with giant molecular clouds. 



ing or even turning-over of the MF. This effect is most pro- 
nounced for initially mass-segregated clusters, and would be 
amplified by the later dynamical cluster evolution, as pre- 
sented in BM03. Although their results can not be straight- 
forwardly combined with the BM03 results (due to differ- 
ences in model setups), both studies suggest even further 
enhan cement of the effects studied in this paper. 

In ILamers et al.l ([2006,1 we constructed simplified evo- 
lutionary synthesis models for solar metallicity, based on 
the GALEV models and the results from BM03. The main 
simplification made concerned the description of the chang- 
ing (logarithmic) MF, which we modelled with fixed slopes, 
but a time-dependent lower mass limit (i.e. assuming that 
only the lowest-mass stars are removed from the cluster, 
while higher-mass stars might only be removed by stellar 
evolution). We scaled our models to match the total mass 
in stars with M < 2 Mq with the BM03 simulations. 

Recently, this appro ach has been improved by 
iKruiissen fc Lameri (|2008D by incorporating the effects of 
stellar remnants for clusters of different initial masses and 
different total disruption times for a range of metallicities. 
They showed that the presence of stellar remnants plays 
a dominant role in the mass evolution of the clusters and 
therefore also in the evolution of the mass-to-light ratio. 
They also found that metallicity affects the colour evolution 
of the clusters, not only by the difference in the colours of 
the stars, but also by influencing the cluster dynamics due 
to the sensitivity of stellar mass and remnant formation on 
metallicity. They determined c olours and mass- to-light ra- 
tios for a range of metallicities. iKruiissenI (|200^ compared 
these predicted mass-to-light ratios with the observed ones 
for cluster samples in different galaxies (Milky Way, Cen 
A, M31 and LMC) and found that the effects of mass seg- 
regation (and the associated preferential loss of low-mass 
stars) can explain the observed range much better than the 
range predicted by s t andar d SSP models. As the models of 
IKruiissen fc Lameri (j2008[ ) are based on the simplified as- 
sumption of only the lowest-mass stars being removed from 
the cluster, they can be improved by models in which the 
mass function changes in a physically more realistic man- 
ner, i.e. by changing the slope of the (logarithmic) mass 
function as derived from dynamical N-body simulations. 
This is the purpose of this paper. 

We describe our input physics in Sect. |2l in particular 
we reanalyse the data presented in BM03 to derive formu- 
lae parametrising the changing mass function (Sect. 12. 3p . 
In Sect. [3] we present our new evolutionary synthesis mod- 
els, and discuss the implications they have for mass-to-light 
(M/L) ratios and cluster age determinations from observa- 
tions. I n Sect, m we present a comparison with prev ious 
models (ILamers et al.ll2006l : lKruiissen fc Lamersl[2008l ) and 
investigate the impact of model uncertainties (fit uncertain- 
ties, initial-final mass relations, isochrones). We finish with 
our conclusions in Sect.|5l 



2. Input physics 

2.1. N-Body simulations by BM03 

BM03 carried out a parameter study for dynamical evo- 
lution of clusters dissolving in a tidal field. They studied 
clusters with a range of particle numbers (8k - 128k, i.e. 
a range in cluster mass) on circular and elliptical orbits at 
different Galactocentric distances (i.e. strengths of the sur- 
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rounding gravitational field). They accounted for mass lost 
due t o stellar evolution (using fit formulae bv'Hurle v et al.l 
|2000|) . two-body relaxation and the external tidal field. 

Th ey initialised their clusters with a universal iKroupal 
()200H ) IMF, which is of the form: 



^(m)dr 



^■^dm 



m < O.5M0 
TO ^ 0.5Mq 



(1) 



with masses in the range O.IM0 ^ to ^ ISMq. This rather 
low upper mass limit was chosen to account for the uncer- 
tain kick velocities of neutron stars (or equivalently, their 
ejection probability from the cluster, see BM03 for details) . 

The MFs provided by BM03 are for single stars, where 
dynamically created binaries are resolved in their compo- 
nents. They give the MF for the whole cluster (i.e. for all 
stars within the tidal radius) . This MF compares well with 
the MF around the half-mass radius of the cluster, as shown 
by BM03. 

BM03 do not take into account primordial mass segre- 
gation and primordial binaries. However, primordial mass 
segregation is found to even furt her increase the change s 
in the MFs found by BM03 fsee iBaumgardt et al.|[2008[) . 
Primordial binaries seem to have little impact on the star s 
evaporating slowly from a cluster (see lKiipper et al. | l2008( ). 
but enhance the number of stars violently ejected during 
strong binary interactions. However, the latter are still only 
a small fraction of the stars leaving the cluster, hence we 
expect little changes in our conclusions if simulations with 
primordial binaries would have been included in our stud- 
ies. 

BM03 do not include a possibly present int ermediate- 
mass b lack hole (IMBH) in the cluster. Recentlv. iGill et all 
(|2008( ) found that the presence of an IMBH reduces mass 
segregation in the centre, which might also influence the 
mass loss from star clusters, although this has still not been 
shown. In addition, the exist ence of IMBHs in star cluster s 
is still under debate (see e.g. iMaccarone fc Servillatl[200il . 

2.2. Total cluster disruption time and the total cluster mass 

Similar to BM03 we will identify the "total cluster disrup- 
tion time" with the time when only 5% of the initial cluster 
mass remains bound. In order to avoid confusion, we will 
specifically label this time tg^%, i.e. the time when the clus- 
ter has lost 95% of its initial mass. However, we will provide 
our models for ages up to the point where a cluster with ini- 
tially 10^ M0 has lost all but 10^ M0 of its luminous mass 
(or to a maximum age of 16 Gyrs, whichever occurs first). 
This termination age of the cluster models is ~ 20 - 26% 
longer than the cluster disruption time tg5% (for models 
with termination ages < 16 Gyr, see Fig. [51 bo t tom p anel) . 

As pioneered by IBoutloukos fc LamersI (|2003t ) and 
iLamers et all (l2005bh . we will use t4 = t^°*'''(M = IO^Mq), 
the total disruption time of a 10'* Mq star cluster, as a 
rough proxy to characterise the strength of the gravita- 
tional field surrounding a cluster: the stronger the field the 
faster the cluster will dissolve, and the shorter t^. Using the 
implicit equation 



j.total 



(Mi) = ti 



IO^Mq 



Mcv (^4 



(2) 



for the total disruption time tjg'*'' (see lLamers et al.ll2005al ). 
t4 can be translated into the total disruption time of clus- 
ters with arbitrary initial mass Mj. For example, a gravita- 
tio nal field chara,cterise d by t4 — 1.3 Gyr (the value found 
bv lLamers et al.ll2005bl for the Solar Neighbourhood) leads 
to complete disruption of a cluster with lO'^ M0 within 
approx. 300 Myr, while a 10^ Mq cluster would survive 
for 22.6 Gyr. /icv(i) describes the fraction of the mass 
that the cluster would have at time t if stellar evolution 
was the only mass loss mech anism, and 7=0.62 (as de- 
termined froin obser vations bv IBoutloukos fc Lamersl[2003l : 
iLamers et al.|[2005bl . a nd in agreement with N-bo dy simu- 
lations, see BM03 and lGieles fc Baumgardl]|2008f) . 

The total cluster mass as a functi on of the fractfonal age 
t/t95% is derived from formula (6) in ILamers et al.l ((200l3) 
(who also show the good agreement with the data from 
BM03): 



t 



95% 



1/7 

[Mcv(Wor-0.05''] ^ (3) 



with Mi = 10^ Mq the initial cluster mass. The stel- 
lar evolution part of this equation was taken directly from 
the GALEV models used in the remainder of this work (for 
details see below). 

Since the total disruption time of clusters in a given 
environment (e.g. tidal field) depends on the initial cluster 
mass, Eq. |3]can also be used to calculate the initial cluster 
mass for an observed present-day total mass and adopted 

t95%- 

The mass fraction in stellar remnants is taken from 
BM03 (then formula (16)): 



/rem(i)=/.^m(i)+0.18. 



t 



^95% 



0.16 



t 



^95% 



(4) 



with f^°^ the mass fraction in stellar remnants from stel- 
lar evolution only (i.e. without dynamical cluster evolution 
effects) taken from our GALEV models, and the other two 
terms describe the increase in the mass fraction of the rem- 
nants due to the preferential loss of low-mass non-remnant 
stars. 

The luminous mass is then: 



Mlumit) - Mtot(t) ■ (1 - /rom(t)). 



2.3. Parametrising the changing mass function 



(5) 



Throughout the paper we look at the logarithm of the 
logarithmically binned mass function (MF). Hence, for a 
Sal peterl (|1955D MF, the power-law index -2.35 becomes a 
linear slope of -1.35. 

In order to parametrise the changes in the (logarithmic) 
mass function, we 

— took the MF data from BMG3 

— divided them by the IMF (by doi ng so we rem ove the 
power law break at 0.5 Mq of the lKroupall200ll IMF) 

— skipped the 2 highest not-empty mass bins (as those are 
affected and partially emptied by stellar evolution), and 

— fitted the remainder with a piecewise power law, inde- 
pendently for every simulation and age. 
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We tried fitting slopes and break points simultaneously, 
finding the results for the break points scattering in the 
range 0.25 - 0.35 M0. This scatter was found to be uncorre- 
lated with any other quantity, confirming it to be of random 
nature. Hence, we chose a double power-law with a break 
point fixed to 0.3 Mq, and only fitted the slopes below and 
above this break point independently for every simulation 
and age. In Fig. [1] we show the slopes of the changing MF, 
relative to the slopes of the IMF (i.e. a{t) — a(0)). The 
points are derived by fitting the data from BM03, using 
the abovementioned scheme, and then the median in bins 
of A{t/tg^%) = 0.025 was taken. The error bars represent 
the 16% and 84% percentiles of individual data points in 
each bin, equivalent to la ranges for a Gaussian distribu- 
tion. These data are grouped according to their disruption 
time t95%: tg^% ^ 3.5 Gyr = filled orange diamonds, 3.5 
Gyr < tg5% ^ 6 Gyr = open blue triangles, 6 Gyr < tg5% 
^ 10 Gyr = filled green circles, and 10 Gyr < tg5% = open 
black squares. 

This fit formula, in conjunction with thelKroupa! (j200l[ ) 
IMF, leads to a 3-coniponent power-law MF with break 
points at 0.3 Mq and 0.5 Mq and time-dependent slopes. 

The time evolution of the changing MF slopes can be 
expressed as: 

a{t) = q;imf(0) 

+ai ■ X + a2 ■ x"^ + as ■ + a4 ■ x'^ (6) 
+b-x-t 

where "t" is the cluster age in Myrs, while "x" = t/tg5% 
is the fractional cluster age in units of its total disruption 
time. The best-fit coefficients are provided in Table[TJ These 
fit parameters have a very high formal accuracy, due to 
the large number of data points used for the fit. However, 
the spread of N-body models around our best fit is the 
dominant source of uncertainty (see below and Sect. 14.1]) . 
We therefore omit the formal fit uncertainties in Table [T] 

In Fig. [1] we overplot our fit formulae for 4 disruption 
times t95% (1 Gyr, 5 Gyr, 8 Gyr and 30 Gyr), representative 
for the chosen grouping in disruption time. 

We restrict the fitting to ages ^ t95%, as in many cases 
clusters with larger ages do not contain enough stars to de- 
termine the MF slopes with reasonable accuracy. However, 
the general trends continue beyond tg5% , following the fit- 
ted relation further on, allowing an extrapolation for ages 
> tg5% (see Fig.H]). 

In addition, we take only simulations into account which 
started with 32k or more particles (these simulations have 
total disruption times in the range 2.3 ~ 25.5 Gyr), as many 
simulations with lower particle numbers show substantial 
uncertainties in the determined MF slopes. 

We would like to emphasize that considering all sim- 
ulations with 16k or more particles or with 64k or more 
particles yield fitted slopes which deviate from the 32k re- 
sults only by less than ± 0.1 for ages up to at least 1.3 x 
^95%- Considering also 8k simulations or only the 128k sim- 
ulations yields larger deviations, due to large run-to-run 
scatter and small number of data/coverage of parameter 
space, respectively. 

On average, the spread of the BM03 simulation results 
around the fitted relation Eq. [S] is of the order of 15%, 
as shown in Fig. [5] and [31 For ages ^ 1/3 x 195% these 
relative deviations are larger, however, at those times the 
absolute deviations of the data from the fit formulae are 
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Fig. 1. Dependence of MF slopes (top panel: for stellar 
masses ^ 0.3 M0, bottom panel: for stellar masses > 
0.3 Mq) on the fractional age t/tg5%. Symbols represent 
median slopes (and 16% and 84% percentiles uncertainty 
ranges) derived from individual data from BM03 runs, and 
binned in intervals of A{t/tg-^%) = 0.025. The data are 
grouped according to their disruption times: tg^o/^ ^ 3.5 Gyr 
(filled orange diamonds), 3.5 Gyr < tg5% ^ 6 Gyr (open 
blue triangles), 6 Gyr < tg5% ^ 10 Gyr (filled green cir- 
cles), 10 Gyr < tg5% (open black squares). Smooth fines 
represent the fit formula ([6]) for the respective age ranges: 
^95% = 1 Gyr (orange dot-dot-dashed fine), 5 Gyr (blue, 
dotted), 8 Gyr (green, dashed) and 30 Gyr (black, sofid). 
Shown is the difference between the time-dependent slope 
of the MF and the slope of the IMF (i.e. a{t) - a(0)). 



small (of the order of Aslope = 0.02-0.03). The impact of 
this uncertainty will be further discussed in Sect. [4TT] 

Although the BM03 simulations are performed for 
a metallicity Z=0.001 (using the fitting formulae from 
iHurlev et al] l2000tl . we will use Eq. [H] for all metallici- 
ties, assuming the metallicity to - at most - introduce 
se cond-order effects on cluster dynamics. This is supported 
bv lHurlev et al.l (|2004l ). who find metallicity effects largely 
cancelling each other, resulting in a weak overall metallic- 
ity dependence of cluster dynamics (although details are 
metallicity-sensitive) . 

2.4. The GALEV models 

The GALEV model s a re extensively described in 
ISchulz et all (|2002[) . lAnders fc Fritze-v. AlvenslebenI 
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Fig. 2. Dependence of MF slope on the age (in Myr) of the cluster for stellar masses ^ 0.3 M0 and for 4 different 
fractional ages t/tg5%: t/tg5% = 0.15-0.25 (top left panel), t/tg5% = 0.45-0.55 (top right), t/t95% — 0.75-0.85 (bottom 
left), and t/t95% — 0.95-1.05 (bottom right). In each panel, the solid line represents the fit formula ^ for the appropriate 
t/t95%, the dashed lines represent the median spread ranges (85% and 115% of fit line). Points are the same as Fig. [TJ 
Please beware of the very different vertical axes. The y-axes are chosen to go from to 2xmaximum fit value. 
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Fig. 3. Same as Fig. [51 but for stellar masses > 0.3 Mq. 
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Table 1. The best- fit coefficients for Eq. ©. The middle 
column gives the coefficients for the low-mass end of the 
MF, for masses ^ 0.3 Mq. The right column gives the co- 
efficients for masses > 0.3 Mq. 



coefficient 


m sC 0.3 M0 


m > 0.3 Mq 


ai 


-0.1345 


0.08389 


a2 


1.7986 


1.9324 


as 


-1.8121 


-0.4435 




1.2181 


0.734 


b 


3.215e-6 


-1.4143e-5 



(|2003l ) and iBicker et al.l (|2004l ). Here we give only a brief 
summary of the relevant input physics used. 

The GALEV models used in this work are based on 
isochrones fr o m th e Padova group, first presented in 
iBertelli et al.l ()1994D . and subsequently updated to in- 
clude the TP-AGB phas43. This update, although not 
documented in a refereed publication, was made publicly 
available approx. 199 9, and treats t he T P-AGB phase 
as later described in iGirardi et al.l (|200G[ ). As we will 
mainly concentrate on the evolution of old stellar clus- 
ters, the Padova isochrones are preferred over the Geneva 
isochrones (Schaller et al. 199§). We want to emphasise 
(and stimulate the various groups of stellar evolution mod- 
ellers) that th e isoch rone sets by IBertelli et al.l (| 19941) and 
ISchaller et al.1 (|1992D (and associated papers) are the only 
available isochrones which cover stellar evolution (in a con- 
sistent way) up to its final stages as well as a mass range 
up to ~ 120 M© required to correctly model ongoing star 
formation in galaxies. For further discussion of this point 
see Sect. 1131 

For consistency with the BM03 simulations we use a 
iKroupal (2001) IMF. 

At each age, the time-dependent MF is evaluated from 
equation ^ for the requested total disruption time tg5%. 
To each star from the iso chrones an appropriate sp ectrum 
from the BaSeL library (jLeieune et al.lll997l Il998t ) and a 
weight according to the time-dependent MF are assigned. 
Then the integrated spectra is obtained by summing up the 
contributions from the individual stars. Here, we assume 
a well-populated MF, hence any stochastic effects due to 
small number statistics, especially at the high-mass end of 
the MF, are neglected. Hence, we model an average star 
cluster. 

The treatment of stochastic effects would go well be- 
yond the scope of this paper. The impact of such stochastic 
effects was studied in depth by Cervifio and collaborators 
(see e .g. 'Cervifio & MoM"200Z; 'Cervifi o fc Luridianal[200l 
|2006| ) and Fagfollni et al. (2007), who found these effects to 
be strongly age- and wavelength-dependent. The strongest 
impact is found for red passbands, which are dominated by 
few red supergiants (young clusters) or very bright upper 
RGB and AGB stars (intermediate-age clusters). The ef- 
fects become smaller for older ages, for non-dissolving clus- 
ters. However, the decreasing number of stars with age in 
our dissolving cluster models likely counteracts this effects. 
We therefore caution users about applying our models to 
a single cluster. The models represent average star clusters 



^ The models of the Padova group are available at their web- 
page: http:/ /pleiadi. pd.astro.it/ 



with the given parameters, hence need to be applied to a 
whole star cluster system. 

More generally, small number statistics is the likely ori- 
gin of the scatter seen in Figs.|2|and|21 However, as we use 
19 BM03 models with a variety of parameters (i.e. total 
masses and dissolution times) to model the dissolution we 
can describe the average cluster dissolution. The impact of 
the spread seen in Figs. ID and |3| will be discussed in more 
detail in Sect. |4?T1 

Due to computational restrictions we only calculated 
individual models for MF slopes with 2 decimal places. If 
at any given age the MF slopes were identical to within 
these 2 decimal places with the MF slopes of a previously 
computed model we reused this older model. Due to this 
finite step-size, some cluster colours exhibit small jumps of 
the order of ~0.001 mag (up to 0.004 mag for the most 
extreme cases). 

The spectrophotometry is normalised to a luminous 
cluster mass as described in Sect. 12.21 

We calculate models for t95%: 

- in the range of 100 - 900 Myr: in 50 Myr steps 

- in the range of 1 - 16 Gyr: in 500 Myr steps 

- for t95% = 18, 20, 25, 30, 40, 60, 100, 150 and 200 Gyr 

and for metallicities (limited by the metallicities provided 
by the Padova isochrones): 

- Z=0.0004 <-> [Fe/H] = - 1.7 

- Z=0.004 ^ [Fc/H] = - 0.7 

- Z=0.008 ^ [Fc/H] = - 0.4 

- Z=O.O2==Z0 ^ [Fe/H] = 0.0 

- Z=0.05 [Fe/H] = + 0.4 

For tg^% > 200 Gyr, wi thin a Hubble t ime the MF slopes 
deviate from the universal iKroupal (|200lh IMF by less than 

0. 005. For conditions similar to the Solar Neighbourhood , 

1. e. t4 — 1.3 Gyr as determined bv iLamers et al.l l2005bl 
the range in total disruption times corresponds to a cluster 
mass range 160 - 3.4 X 10'^ Mq (i.e. a 160 Mq cluster needs 
100 Myr to totally disrupt, a 10^ Mq cluster needs 1.3 Gyr, 
and a 3.4 x 10^ Mq cluster needs 200 Gyr). For the SMC, 
with t4 — ^10 Gyr, the range in total disruption times 
would correspond to a cluster mass range 6 - 1.25 x 10^ 

Mq. 

The data are made publicly available at our web- 
pages http://www.phys. uu.nl/'^anders/data/SSP_varMF 
and http://data.galev.org. We provide the user with inte- 
grated cluster magnitudes in a variety of passbands plus 
cluster masses (total mass, luminous mass and mass in stel- 
lar remnants) for each of the models. Integrated spectra are 
available upon request. 
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3. Results and implications 

In this section we present our new models for solar metal- 
licity (unless otherwise noted) and discuss its implica- 
tions. Models for all metallicities are available at our web- 
pages http://www.phys. uu.nl/^anders/data/SSP_varMF 
and http://data.galev.org. We would like to emphasize that 
the absolute values of our models (and therefore also the 
results and implications discussed in this section) depend 
on our choice of isochrones and other input physics. In Sect. 
IHwe discuss some of these uncertainties. The main results, 
the systematic differences induced by the preferential mass 
loss, are hardly changed. 

3.1. Photometry 

The photometry for the new models is shown in Fig. 3] 
The colours are shown as differences between the new mod- 
els with changing MF and the standard models with fixed 
(initial) MF. For illustrative purposes, in the bottom left 
panel the absolute values for the V-I colour are presented. 
The V-band magnitude evolution (bottom right panel) is 
given in absolute values for a 10^ Mq cluster. 

The V-band magnitude evolution shows the stellar evo- 
lution fading line as bright limit for the new models, which 
they follow for young ages, when the effect of mass loss is 
not yet pronounced. Already after ~10% of their respec- 
tive total disruption times, the new models have evolved 
0.1 mag away from the fiducial fading line, due to the loss 
of stars. At ^80% (± 10%, depending on the model) of 
their respective total disruption times, the new models are 
1 mag fainter than standard models predict, due to the loss 
of stars. 

Except for the very earliest stages of cluster evolution 
(first few Myrs), the flux in passbands redder than the V 
band is dominated by stars initially more massive than the 
main contributors to the flux in bluer passbands. This is 
caused by the flux in the red passbands being dominated 
by red (super)giants, which are more luminous than the 
stars at the low- mass end of the main sequence (MS), even 
after taking the larger number of low-mass MS stars due to 
the IMF into account. For a changing MF due to dynamical 
evolution, the contribution from low-mass MS stars is even 
further reduced. 

The dominant source of flux contribution in passbands 
including and bluewards of the V band is a strong function 
of time: in early stages, the flux is dominated by mid-MS 
stars (the evolution through the Hertzsprung gap is too 
fast to significantly contribute). As the cluster ages, the 
MS turn-off (MSTO) shifts successively redwards through 
the filters, ever increasing its contribution to the band's 
flux. However, the relative contribution of mid-MS stars 
and MSTO stars is also strongly dependent on the MF, 
hence is dependent on the total disruption time of our mod- 
els. 

As the selective mass loss preferentially removes the 
least massive stars from the cluster (and therefore its inte- 
grated photometry), it causes the cluster to generally be- 
come redder than the standard models without cluster dis- 
solution (i.e. with infinite total cluster disruption time). 
The MF evolution and the resulting reddening speeds up 
while the cluster approaches its final disruption, leading to 
the steep colour evolutions towards the end of a cluster's 
lifetime, as seen in Fig. 21 



Two exceptions have to be noted: 

— the colour U-B (and similar colours) becomes bluer than 
the standard models for total disruption times shorter 
than ^ IGyr. At these ages the B band is entirely dom- 
inated by mid-MS stars, while bluer bands contain con- 
tributions from the higher MS stars and the MSTO 
stars. As the mid-MS is stronger depopulated than the 
upper MS and MSTO due to the dynamical cluster evo- 
lution, the colours become bluer. For longer total disrup- 
tion times the mid-MS is not sufficiently depopulated 
to notice the effect until the B band gets contributions 
from the MSTO stars. Redder passbands do not show 
this effect as they contain contributions from bright red 
(super-)giants. 

— for colours like V-R and V-I and ages 6 Gyr the mod- 
els get slightly bluer than the standard models for total 
disruption times <; 10 Gyr. Likely, this comes from the 
strong depopulation of the lower MS (and the standard 
IMF containing a large number of stars at low masses), 
which leaves an imprint even though a single lower-MS 
star is 3-4 mag fainter than an RGB/AGB star of simi- 
lar temperature. Redder passbands do not show this ef- 
fect as the magnitude difference between lower-MS stars 
and RGB/AGB stars increases with increasing wave- 
length/decreasing temperature, and the total contribu- 
tion from MS stars decreases. 

For long total disruption times, the maximum colour 
deviation of our models for dissolving clusters from the 
standard models decreases with increasing total disruption 
time. This comes from several effects at once: 

— for such old ages, the MF covers only a narrow mass 
range in both cases 

— the integrated cluster flux is dominated by the upper- 
RGB/ AGB stars, as they are significantly brighter than 
the MSTO region (the magnitude difference between 
upper- RGB/ AGB and MSTO increases with time), re- 
sulting in an even narrower "effectively visible MF" 
range 

— the temperature range these stars cover is significantly 
smaller than at younger ages, resulting in a lower sen- 
sitivity of the colours to the exact distributions of stars 
along the isochrone. 

However, the changes in mass/absolute magnitude (see 
Fig.m bottom right panel) and M/L ratio (see next section) 
are significant in all cases. 

The increasing maximum colour deviation of our dis- 
solving clusters models from the standard models for blue 
passbands and increasing total disruption times < 1 Gyr 
originates from the redward shifting of the MSTO through 
the filters. 

3.2. Mass-to-light ratios 

The low-mass stars preferentially removed in the course of 
cluster dissolution have mass-to-light (M/L) ratios which 
are higher than the M/L ratio of the average cluster star. 
On the other hand, as shown by BM03, the fraction of 
(non-luminous) stellar remnants in dissolving clusters is en- 
hanced w.r.t. the standard models. These effects partially 




Fig. 4. Solar metallicity models with the changing mass function treatment, following equation ([6]) with t95% in the range 
100 Myr - 100 Gyr (for clarity only models at 100, 200, 400, 800 Myr and 1, 2, 4, 8, 10, 12, 14, 16, 18, 20, 40, 100 Gyr 
are shown), going from left to right, respectively. Upper left panel: U-B colour, top right: B-V, middle left: V-I, middle 
right: V-K. For these colours, the differences between the new models and the standard models are displayed. Bottom 
left: the absolute values for V-I, bottom right: V-band magnitude evolution of a 10^ Mq cluster. Thick lines represent 
the standard models. 



cancel each other and lead to the time-dependent M/L ra- 
tios shown in Fig.[5l It demonstrates that for each total dis- 
ruption time the M/L ratio of our dissolving cluster models 
is systematically lower than a standard model would sug- 
gest for the majority of a cluster's total disruption time. For 
the final up to '^16% of a cluster's total disruption time, 
the M/L ratio can become enhanced w.r.t. the standard 
models (see Fig. [51 bottom panel). This is caused by the 
increasing fraction of stellar remnants inside the cluster. 

In Fig.[6]we present the dependence of the V-band M/L 
ratio on the present cluster mass. The top panel shows how 
this relation evolves with cluster age at a field strength (i.e. 
location in a galaxy, characterised by t4, the total disrup- 
tion time of a 10''' Mq cluster, as described in Sect. 12. 2p 



representative for th e Solar Neighbourhood, as found by 
iLamers et al.l (|2005a|) . As the clusters evolve, generally the 
M/L ratio increases due to stellar evolution. In addition, 
the lowest-mass clusters eventually disrupt (and drop out 
of this plot). The highest-mass clusters lose mass, but still 
have M/L ratios close to the canonical value from stellar 
evolution. Intermediate-mass clusters suffer from the im- 
pact of the changing mass function which reduces their M/L 
ratio significantly compared to the canonical value. A few 
cases of enhanced M/L ratios can be seen for clusters close 
to final disruption (at the low-mass end of the curves). 

In Fig. [S] (top panel) we overplotted data of 
young (ages < 1 Gyr) LMC a nd SMC clusters by 
iMcLaughhn fc van der MareH (pOOSh (labelled "MM05" , for 
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details of this dataset see following subsection) and Larsen 
and collaborators (labelled "L arsen" ) in 4 spi ral and irreg- 
ular g alaxies, (see [Larsen fc R ichtlcr 200^and lLarsen et al.l 
|2004|) . Out of the 13 clusters in these samples, 8 have M/L 
ratios consistent with our models for their respective ages 
(within their la uncertainty ranges). The remaining 5 clus- 
ters all have too high M/L ratios for their respective ages. 
2 clusters are very young (^^10 Myr), hence could be out 
of equilibrium during their gas expulsion and readjustment 
phase, and their veloci ty dispersions might not t race their 
dynamical masses fsee lGoodwin fc Bastiai3l2006l) . The de- 
viations of the remaining clusters might be indications of 
errors in the models, or could be signs that the age determi- 
nation is uncertain or the velocity dispersions are seriously 
affected by the orbital motions of binaries (or other sys- 
tematic observational effects, like macroturbulence in the 
stellar atmospheres or instrumental resolution). 

The middle panel shows the M/L ratio in the V band 
as a function of the present-day mass of 12 Gyr old clus- 
ters, for a range of gravitational field strengths (i.e. typical 
disruption times 14). Within each line, the cluster's total 
disruption time goes from 10 Gyr (low mass end; clusters 
with shorter total disruption times have been disrupted by 
an age of 12 Gyr) to 200 Gyr (upper end of available total 
disruption time range). For example: a cluster located at a 
position in a galaxy characterised by a field strength with 
t4 = 1.3 Gyr (i.e. the blue line, corresponding to the en- 
vironment in the Solar Neighbourhood), observed now (i.e. 
at an age of 12 Gyr) with a mass = 10^ Mq is expected to 
have a M/Ly of ^^4, while a cluster with a mass — 10"* M0 
has a modelled M/Ly of ~2.3. 

The bottom panel shows the impact of metallicity on 
the M/L ratios: with increasing metallicity the M/Ly ra- 
tio based on stellar evolution increases. Therefore all curves 
are stretched to reach higher M/Ly ratios for higher metal- 
licities, while the shape of the curves is largely unaffected. 
At the high-mass end, all curves level off to their respective 
values determined by stellar evolution alone. 

By comparing the top and bottom panels of Fig. [6l both 
for our new models as for the non-dissolving standard mod- 
els the well-known age -metallicity degeneracy is apparent 
(see e.g. iWorthevlfTOOl . 

As the GALEV code (like most other evolutionary syn- 
thesis codes) is not capable to directly deal with stochastic 
effects (especially not with the stochastic effects inherent 
to the selective mass loss caused by dissolution), we em- 
ploy the following estimation scheme: 

— We assume, that the majority of stochasticity originates 
from evolved stars (mainly RGB, AGB stars). 

— As a function of age, we determine the relative number 
of evolved to unevolved stars (i.e. MS stars). 

— From this relative number we determine the total num- 
ber of evolved stars for clusters of different masses, and 
the stochastic scatter (i.e. the square root of the total 
number of evolved stars). 

— We determine average properties of the evolved stars 
(mean effective temperature, mean log(g) and mean lu- 
minosity) . 

— We multiply the stochastic scatter with the spectrum of 
the mean evolved star, and add resp. subtract this from 
our standard spectrum. 

From this approach, we estimate the effect of 
IMF stochasticity on the M/L ratios to be roughly: 



15/5/1.5 % uncertainty for clusters with total mass 
10,000/100,000/1,000,000 Mq, respectively. This test was 
only done for the standard, not depopulated IMF. The ef- 
fects for the depopulated MF will be smaller, as for the 
same total mass, the number of giant stars will be larger. 

Our results show fe a tures similar to those presented 
by 'Kruiissen fc LamersI ((2008') and applied by 'KruiissenI 
((2008,). However, systematic differences are present, inher- 
ent to the underlying assumptions, and discussed in Sect. 

3.2.1. Comparison with observations 

In the following section we want to compare our new models 
with old globular clusters (and other old massive stellar 
systems) in the Milky Way and other galaxies. This is a 
first step to validate our models. 

Fig. [7] compar es our new models with obs erva- 
tional data from .McLaughlin fc van der Mareil (|2005D and 
iMieske et all (|2008D . The colour coding of the data refers 
to their metallicity: 

- blue = [Fe/H] < -1.2 = "MP" (metal-poor) 

- green = -1.2 < [Fe/H] < -0.55 = "IM" (intermediate 
metallicity) 

- red = -0.55 ^ [Fe/H] < -0.2 = "MR" (metal-rich) 

- magenta = -0.2 [Fe/H] = "SO" (around solar). 

These ranges were chosen to be consistent with the 
metallicities of the Padova isochrones/our models. 3 objects 
from the cited samples are not shown in these plots due 
their high M/L ratios: LMC-NGC2257 an d MW-NGC6535 
from lMcLaughlin fc van der Mareil (|2005D . which have M/L 
ratios of the order of 8-10 with error bars of the order of 
4-5, and the Virgo cluster o bject S999 (with M/L = 10.2, 
from the iMieske et al.l |2008| sample) , which might have a 
genuinely high M/L ratio. 

Overplotted are 2 bundles of models for metallicities 
in the range from [Fe/H]= - 1.7 to 0.0 (as is appropriate 
for the shown observational data), for cluster ages of 12 
Gyr, and local tidal field strengths with 14 = 5 Gyr (up- 
per/left branches of models of a given metallicity, repre- 
sentative for halo clusters) and = 300 Myr (lower/right 
branches, representative for strong dissolution). The main 
point of the comparison is to show the range of M/Ly val- 
ues reachable with our models. As can be seen in Fig. [71 
most Galactic GCs have M/L values compatible with our 
predictions and many, especially low-mass ones are below 
those of the standard isochrones. The estimated uncertain- 
ties of few per cent, as estimated above for clusters in this 
mass range, are not sufficient to bring the observations into 
agreement with the standard predictions from stellar evolu- 
tion alone. We take this as clear evidence for cluster evolu- 
tion/dissolution and that our evolving cluster models are a 
clear improvement over standard isochrone fitting for GCs. 

Data for the Milky Way a nd the LMC are taken from 
[McLaughlin fc van der Mareil (|2005t ). the most extensive 
homogenised compilation of star cluster M/L ratios (and 
providing also other star cluster properties) for these galax- 
ies. The m ajority of the data for the Milky Way originates 
from lPrvo r fc Meylan 1993. P rvor fc Mevla n (1993) find a 
weak correlation of M/L ratio with mass, consistent with 
our models, but with large scatter and uncertainties (on 
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average about 50-60%) and M/L ratios outside the acces- 
sible range of our models for some of their sample clus- 
ters (fMcLaughlin & van der Marel 200J^ do not elaborate 
on this dependence) . They find no significant correlation of 
M/L ratio with the distance of the cluster from the Galactic 
Center or the Galactic plane, contrary to what might be 
expected from the BM03 simulations and our models (al- 
though the mixture of clusters with different masses at dif- 
ferent Galactocentric radii, i.e. experiencing different tidal 
field strengths, could erase any such signature). However, 
the present-day cluster position within the Galaxy is likely 
less important for the total disruption time (and there- 
fore the M/L ratio evolution) than the perigalactic distance 
and the number of past disk passages, which are unknown 
for most clusters. In addition, stochastic effects of the MF 
could induce additional scatter. The error bars are too large 
to find a clear trend of M/L ratio with metallicity. 

Of the 52 old clusters all but 8 are consistent within 
their la ranges with models for [Fe/H]= - 1.7 or [Fe/H] = 
- 0.7. All of these 8 clusters have significantly too low M/L 
ratios. NGC 241^ and NGC 4590 both have metallicities0 
below [Fe/H]= -1.7, the lowest metallicity for which we can 
provide models. Those clusters could possibly be explained 
by models of even lower metallicity. For the other clusters 
(NGC 5272, NGC 5286, NGC 5904, NGC 6366, NGC 6715, 
NGC 7089), no immediate explanation (apart from under- 
estimated observational uncertainties or the impact of the 
unknown perigalactic distance and past disk passages) is 
apparent. However, our new models represent a significant 
improvement: while 6 clusters are not consistent with our 
new models, 21 clusters are not consistent with the stan- 
dard constant-M/L models. 

Another way to analyse the properties of Milky 
Way globular clusters is Fig. [8l where we com- 
pare their dereddened V-I colours (taken from the 
Harris catalogue) with their M/L y ratios (as given by 
iMcLaughlin &: van der Marelll2005[ ). We over plot our mod- 
els for a cluster age of 12 Gyr. The models with the 
longest disruption time are equivalent to the standard/non- 
dissolving model (marked with the red asterix). For de- 
creasing disruption time, the models' M/L ratios drop, be- 
fore drastically increasing again at the final stages of dis- 
solution. We restrict this analysis to metal-poor clusters 
(-2 < [Fe/H] < -1.2), as the number of higher- metallicity 
clusters with the required data is too low to draw strong 
conclusions. We find good agreement between the obser- 
vational data and our models for these clusters concern- 
ing their M/Ly ratios: the observational data are clearly 
spread out over a wider range than the standard model 
could account for, while our new models cover this range 
much better. However, the observed cluster colours span a 
wide range in V-I (though no colour uncertainties are avail- 
able), which we cannot fully account for with our models. 
Our models might be about 0.05mag too red for the obser- 
vations. This could originate from our choice of isochrones: 
see Sect. 14.31 where we find that other isochrones give re- 



^ A re-analysis of the velocity dispersion of NGC 2419 by 
Baumgardt et al., submitted, indicates that the mass-to- light 
ratio is around 2, which is in good agreement with a canonical 
mass-to-light ratio and no dynamical cluster evolution. 

* Data taken from the Harris catalogue iHarrij (|1996| ). avail- 
able at 

http:/ /physwww. physics. mcmaster.ca/^harris/mwgc.dat 



suits that are bluer than our set of isochrones. However, 
the other isochrones are than bluer than the observations, 
by again ~0.05mag. Other possible causes include: uncer- 
tainties in the reddening estimates, filter curve mismatch, 
etc. 

Data for massive star clusters in NGC 5128 {— Cen 
A) as well as for massive objects (commonly referred to as 
"Ultra-Compact Dwarf galaxies" = U CDs) in the Virgo and 
Fornax galaxy cluster are taken from lMieske et al. (200^ 
These data include earlier observations bv iReikuba et all 
(|2007D for the star clusters in Cen A, a nd observations by 
a variety of authors for the UCDs (see lMieske et al]|2008l 
for details). While 27 of these clusters are not consistent 
with our new models, 47 clusters are not consistent with 
the standard constant-M/L models. Also for this sample, 
the new models are a significant improvement. 

We therefore conclude that the samples are better de- 
scribed by our new models with preferential loss of low- 
mass stars, and we see ongoing cluster dissolution. In fu- 
ture cluster modelling, this effect has clearly to be taken 
into account. 

Nonetheless, the sample of massive Cen A clusters and 
UCDs shows a very clear and strong trend of increas- 
ing M/L ratio with object mass, especially for metal- 
poor/intermediate metallicity objects (metal- rich objects 
are reasonably well covered by our models, except for the 
Virgo cluster UCD S490). This trend can not be repro- 
duced by our models: for masses larger than wlO"^ Mq 
only 3 out of 12 objects are consistent with our models 
within their respective la uncertainties (one further ob- 
ject is marginally consistent). Such massive systems are 
not expected to be mass segregated due to their large re- 
laxation time, let alone close to disruption (which in our 
models is the only possibility to reach M/L ratios higher 
than the predictions from standard stellar evolution). While 
the models do have inherent sources of uncertainties (e.g. 
the assumed initial-final mass relation for remnants, uncer- 
tainties in the underlying stellar isochrones, which will be 
studied in more detail in Sect. [4]), they are unlikely to raise 
the model M/L ratios sufficiently to accommodate a signifi- 
cant fraction of the currently unexplained observations (es- 
pecially without removing the agreement for objects with 
lower M/L ratios). Two possible explanations for the high 
M/L ratios would remain: either a stellar mass function 
significantly deviating from the universal iKroupal (I200lf) 
IMF (see also lDabringhausen et al.ll2008HMieske fc Kroupal 
I2008D . or dark matter fsee iBaumgardt k Mieskd l2008l for 
how dark matter can explain the high M/L ratios of UCDs). 

3.3. Impact on age determination 

Evolutionary synthesis models are regularly used to derive 
the physical parameters of star clusters (and galaxies) from 
observed spectrophotometry. Derived quantities are age, 
mass and metallicity of the star cluster as well as the extinc- 
tion in front of the star cluster (sec e.g., among many others, 
Bicker ct al. 2002; Kassin et al. 2003[jLnders et ah 2004S 
de Griis et al. 2004; Kundu et all 120051: Ide Griis fc AndeTsI 
i2006l: ISmith et al.i i2007l) . Our GALEV models provide a 
model grid of SEDs in age/metallicity/extinction. The 
"AnalySED tool" (which we developed and tested in 
lAnders et al.ll2004af) compares these model SEDs with the 
observed SED of a star cluster using a algorithm, to 
derive the best-fitting parameter combination and their 
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respective uncertainty ranges from integrated multi-band 
cluster photometry. 

Here we employ the "AnalySED tool" to quantify the 
differences between the true ages of dissolving clusters (with 
a time-dependent MF) and the ages derived using the stan- 
dard evolutionary synthesis models (with a MF fixed to the 
IMF slopes) . We take the cluster photometry from the dis- 
solving cluster models (for a number of filter combinations), 
apply Gaussian noise (with a = O.lmag) to the photom- 
etry in the individual passbands, and analyse them using 
the standard, non-dissolving cluster models. The analysis is 
done for fixed solar metallicity and zero extinction (leaving 
these parameters free to vary would lead to even stronger 
deviations from the standard models and larger uncertain- 
ties, as shown in I Anders et al.ir2004a[ ). For each filter com- 
bination, total disruption time and age we generate 1000 
test clusters, derive their physical parameters and deter- 
mine the mean of the derived ages. The results in terms of 
the ratio between the derived mean age and the true cluster 
age are shown in Fig. [21 

In general, for all models the ages get overestimated 
for a significant fraction of the cluster lifetime (for some 
ages and models the ages can also get severely underes- 
timated). This agrees well with the discussion concerning 
the cluster colours in Sect. 13.11 generally, when the clus- 
ter colours become redder than the standard models, the 
ages become overestimated. A direct comparison is not ap- 
propriate, though, as "AnalySED" uses the whole available 
Spectral Energy Distribution (SED, i.e. the dataset con- 
taining all magnitudes in a given set of filters for a given 
cluster) to determine the model with the best-matching pa- 
rameters, hence differences in different filters can either can- 
cel or amplify each other. 

Datasets including the mid-UV (here represented by the 
ACS HRC F220W filter) show only modest deviations from 
the standard models (Fig.[9l top panels). However, 20% de- 
viations are regularly found. Datasets lacking the mid-UV, 
and especially those including near-IR data, are more sen- 
sitive to the changes in the mass functions (Fig. [9l bottom 
panels). For those datasets, deviations of 50% or even a 
factor ~2-3 are found. 
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Fig. 5. Same models as Fig. HI but M/L ratios. Upper 
panel: time evolution of V-band M/L ratio for the new mod- 
els. Middle panel: ratio of new model's V-band M/L ratio 
and V-band M/L ratio from standard models. The non- 
dissolving/standard model is shown as thick line. Bottom 
panel: characteristic times of the models as function of to- 
tal disruption time: black crosses and line = age at which 
the cluster contains only 100 Mq luminous matter (i.e. ter- 
mination age of model); red open circles = age at which 
the M/L ratio evolution crosses the standard model; green 
filled squares — age at which the M/L ratio of the dissolving 
cluster models is minimal w.r.t. the standard models (i.e. 
the dip seen in the middle panel). Only models for which 
the respective age is smaller than the maximum model age 
of 16 Gyr are shown. 
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Fig. 9. Derived ages of dissolving star clusters with tg5% in the range 100 Myr - 60 Gyr (from left to right: 100, 500 
Myr, 1, 5, 10, 20 Gyr), using standard models, for a number of filter combinations. Upper left panel: ACS HRC F220W 
+ standard UBVRIJHK, top right: ACS HRC F220W + standard UBVK, bottom left: standard UBVK, bottom right: 
standard UBVI. Shown is the ratio between the age derived using standard models and the age the dissolving cluster 
has. 
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4. Validation of the models 

In this section we will investigate several uncertainties in- 
herent to our models, as well as comparing our new models 
with models previously released by our group. 

All the following values are maximum differences from 
the standard models of a given parameter in a given time 
interval, unless otherwise noted. In many cases the maxi- 
mum deviations occur in the final stages of cluster dissolu- 
tion, and for the longest total disruption times < maximum 
model age. If we would have chosen a maximum model age 
of 13 Gyr instead of 16 Gyr, the maximum deviations would 
generally be slightly smaller. 

For the different issues discussed in this section, we also 
publish a few test cases on our webpage, illustrating the im- 
pact of different initial-final mass relations, isochrones and 
parametrisations of the mass function evolution on colours 
and masses/ mass-to-light ratios. For the parametrisations 
of the mass function evolution we select a few disruption 
times for presentation on the webpage, while for the initial- 
final mass relations and isochrones we present only data 
without disruption (i.e. pure stellar evolution) to avoid con- 
fusion. 



4.1. Parametrisation of mass function slope evolution 

As discussed earlier, the mass function slope evolution is 
derived from a subset of N-body simulations by BM03. The 
subset was selected to cover the parameter space well, while 
limiting the impact of low- number statistics (see Sect. 12. 3t . 

The fit to the data of the time evolution of the mass 
function slopes has formally a very high accuracy due to the 
large number of data points. However, as shown in Figs. [2] 
and [31 the N-body models show an intrinsic spread around 
the fitted function. We quantified this spread to have a 
median value ^ ±15% for ages ^ 1/3 tg5%. For younger 
ages, this relative spread is larger, however the median of 
the absolute spread is small, ^ ±0.02-0.03 change in the 
slope. 

We test the impact of this spread by calculating mod- 
els for which the time-dependent part of the mass func- 
tion slope is reduced/increased by 15%. We find, as ex- 
pected, the change in the "high-mass slope" (i.e. for masses 
^ 0.3 Mq) to be of primary importance, while the time- 
dependent contribution from stars with masses ^ 0.3 Mq 
changes the photometry only mildly. 

As the mass evolutions of the cluster (total, luminous 
and remnant mass) were derived independent of the mass 
function evolution, the masses are not affected. 

The impact of this uncertainty on the colors is small: the 
changes induced for the models with the shortest disruption 
time tg5% (i.e. 100 Myr) and the colors with the longest 
wavelength coverage (i.e. V-K) reach w O.OTmag at final 
disruption. These changes decrease rapidly with increasing 
disruption time and decreasing wavelength coverage. 

For ages t ^ tg^% the magnitudes change by ~ 0.15 - 0.2 
mag (with the changes slightly larger for the shortest tg5% 
and red passbands). As the mass is unaltered, this directly 
translates into a change in the M/L ratio by 15 ~ 20%. 

For ages t > tg5% both magnitudes and M/L ratios di- 
verge from the models using the best fitting relation for the 
time evolution of the mass function slopes. Models with a 
weaker time evolution are increasingly brighter and have 
consequently lower M/L ratios. 
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Fig. 6. Top panel: time evolution of V-band M/L ratio 
(for solar metallicity) as a function of present-day clus- 
ter mass for a fixed local gravitational field strength char- 
act eristic for the Solar Neighbourhood (i.e. t4=1.3 Gyr, 
see iLamers et al.l l2005al ). Middle panel: V-band M/L ra- 
tio (for solar metallicity) as a function of present-day clus- 
ter mass for a range of local gravitational field strengths 
for clusters observed at an age of 12 Gyr. Bottom panel: 
V-band M/L ratio as a function of present-day cluster 
mass for a range in metallicity for clusters observed at 
an age of 12 Gyr and experiencing a typical disruption 
ti me t4=1.3 Gyr. In the to p pan el, observations from 
Mc Laughlin fc van der Marell (I2005D fMM 05) and from 
iLarsen fc Richtler 2004 : Larsen et al.ll2003 are overplotted, 
for young clusters with ages < 1 Gyr. 
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Fig. 7. Observations of old ( ~12 Gyr) objects from 
McLaughlin &: van der Marell (|2OO50 (MM05) and 
Mieske et alj l2008l (M08). Top panel: objects with 
[Fe/H] < -0.55. Bottom panel: objects with [Fe/H] ^ 
-0.55. Overplotted are our model V-band M/L ratio as a 
function of present-day cluster mass for the appropriate 
metallicities, at 12 Gyr and for typical disruption times 
t4=300 Myr (lower/right branches for each metallicity) 
and t4=5 Gyr (upper/left branches). Horizontal lines at 
the canonical values of M/L from stellar evolution are 
also shown. They represent the limit of infinite disruption 
times. These models are intended to illustrate the range 
covered by our models, and the inconsistency of most 
obsverations with models of infinite disruption time. 



In Fig. [10] we illustrate the effects of enhancing the MF 
evolution by 15%. Diminishing the effects of MF evolu- 
tion by 15% gives quantitatively similar results with the 
changes w.r.t. the standard models going in the opposite 
direction. The V-K colour evolutions are the most extreme 
cases: Effects become smaller for shorter wavelength cover- 
age and longer disruption times. 

In summary, the uncertainty induced by the spread of 
N-body model data around the fitted time evolution of the 
mass function slopes has an impact on the model predic- 
tions. However, for ages t ^ tg5% the induced uncertainties 
are much smaller than the error one makes by not taking 
into account the effect of preferential mass loss and cluster 
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Fig. 8. Observations of old (~12 Gyr) Milky Way glob- 
ular clusters from the Har r is cat alogue (colours) and 
iMcLaughlin fc van der Marell (|2005f ) (M/L ratios). Shown 
is the relation between observed V-I colour (dereddened) 
vs M/Ly ratio for metal-poor clusters. 
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Fig. 10. Integrated V-K colour (left panel) and M/Ly ratio 
(right panel) for a Kroupa IMF, solar metallicity, our stan- 
dard isochrones, a range of cluster disruption times, with 
the effect of MF evolution enhanced by 15% w.r.t. the stan- 
dard models. Shown are the quantities relative to the re- 
spective quantities of our standard models. Diminishing the 
effect of MF evolution by 15% yields quantitatively similar 
results, however, the changes are in the opposite direction. 
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Fig. 11. M/Ly ratio for a Kroupa IMF, solar metallicity, no 
cluster dissolution and the 3 alternative IFMRs discussed in 
the text. Shown are the quantities relative to the respective 
quantities of our standard models. 



dissolution. In addition, we aim at describing the average 
cluster. 



4.2. Initial-final mass relations 

Our model results, especially the M/L ratios discussed in 
Sect. 13.21 depend on the treatment of stellar remnants. 
The remnants mass is calculated from the progenitor star's 
initial mass and the adopted initial-final mass relation 
(IFMR), which accounts for mass loss during the life of 
the progenitor star and due to the "death" of the star and 
the formation of the remnant. 

The IFMR fo r white dwarfs used in t his wo rk is based 
on the work by IWeidemann fc Koesterl (|l983t) (hereafter 
"Weidemann83" ) . As the IFMR is still uncertain, we tested 
our choice by adopting differ ent IFMRs for white dwarfs, 
na melv bvlWeidemannI (|2000f ) (hereafter "WeidemannOO"), 
bv lKahrai et al l (^2008^ (h ereafter "KahraiOS") and the pre- 
scription bv lHurlev etall (|2000D (hereafter "HPTOO"). For 
the latter one we also adopt their I FMR for neutron stars, 
while for all other IFMRs we adopt iNomoto fc Hashimotol 
fl988[ 

Changes discussed below are w.r.t. our standard IFMR 
Weidemann83. 



We find the IFMR to be of minor influence on the results 
for ages t ^ tg5% : the total mass changes by maximum 2 - 
3%, while the luminous mass changes by 4% and 6% (for 
WeidemannOO/KaliraiOB and HPTOO, respectively). This 
translates into magnitude changes of ~0.05mag and ^0.07 
mag (for WeidemannOO/KaliraiOB and HPTOO, respec- 
tively). The associated effect on the M/L ratios is ~4.5% 
and ~6.5% (for WeidemannOO/KaliraiOB and HPTOO, re- 
spectively) . 

For larger ages t ^ tgQ% the results eventually diverge. 
However, only in the last 5% of a cluster's lifetime the total 
mass differs by more than 10%, regardless of the choice of 
IFMR. 

In Fig. [11] we show the impact of the chosen IFMR on 
the M/Ly ratio for infinite disruption time. 



4.3. Isochrones 

Our choice of isochrone s (i.e. isochrones fro m the Pavoda 
group, first presented inl Bertelli et al.|[l994l with later up- 
dates of the Padova group concerning the TP-AGB phase 
= "updated Padova94") is driven by the following points: 

1. For consistency with GALEV models of galaxies we re- 
quire isochrones which cover the full mass range up to 
high masses (ideally up to ^ 120 M©) to properly model 
ongoing star formation in galaxies. 

2. Likewise models covering a wide range in metallicities 
is desired to consistently model old/metal-poor globu- 
lar clusters and young/metal-rich star clusters formed 
in nearby starbursts, as well as to model galaxies consis- 
tently from the onset of star formation to their present 
stage. 

3. The models should cover all relevant evolutionary stages 
of these stars, especially the very luminous phases (for 
our study especially the TP-AGB phase is of prime im- 
portance, but also early stages like supergiants are im- 
portant) . 

We regret that recent high-quality isochrone c alcula- 
tions (e.g. Girardi et al. 2000; Yi et al. 2001; Cariul oet al.l 
120041: iPietrinferni et al.ll2004 1200a iBertem et al.ll2008l) are 
all focussing on "low-mass" sta rs (maximum up to ~ 10 
Mq, though iBertelh et al.|[2008l announce models up to 20 
Mq for the near future) and/or do not fulfil one or more cri- 
teria mentioned above. However, stars more massive than « 
10 Mq contribute significantly to the chemical enrichment 
and the light of young star clusters and most galaxies. 

The only models fulfillin g all mentioned cr iteria are 
models by the Padova group f Bertelli et al.l 'l994 pl us TP- 
AGB updates) and by the Geneva group ( Schaller et al.l 
[l99l IChaTbonnel et al.|[l99l iSchaerer et al.lll99l . 

As the main focus in this paper is on systems older than 
^ 100 Myr we prefer the updated Padova94 isochrones over 
the Geneva isochrones. 

The alternative solution, to combine isochrones from 
different groups/epochs, was rejected as consistency cannot 
be ensured. 

W e tested solar-metal l icity i s ochrones by Cariulo et al.l 
()2004t ): IPietrinferni et al.l (|2004D : IMarigo et al.l (|200B[ ) falso 
known as "Pisa/GIPSY" , "BASTI" and "new Padova", 
respectively) with respect to the updated Padova94 
isochrones we used in this study, and derived star clus- 
ter models for test purposes. While the Pisa isochrones 
are offset from all other isochrones (they are generally sig- 
nificantly hotter, but are based on more limiting input 
physics) , the other isochrones are in overall good agreement 
with the updated Padova94 isochrones. Small differences 
include: 

— For increasing age, the BASTI main-sequence turn-off 
temperature goes from slightly cooler than the updated 
Padova94 to slightly hotter (by a few per cent). This 
results in an increasing deviation of U-/B-band mag- 
nitudes compared to the updated Padova94 isochrones 
by up to 0.5mag. The new Padova isochrones show 
much smaller deviations ^ O.lSmag in these passbands. 
Contrary, for both BASTI and new Padova, colours like 
U-B or B-V deviate for most of the time by ^ O.lmag, 
and for the majority of time by ^ O.OSmag from the 
updated Padova94 models. 
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— Overall, the RGBs and AGBs in the BASTI and new 
Padova isochrones are hotter than in the updated 
Padova94 isochrones. Especially stars with the highest 
luminosities on the RGB/AGB are treated differently. 
For ages younger than ~1 Gyr, the test models deviate 
significantly, both from our standard model as well as 
from each other. For ages J52 Gyr, the BASTI and new 
Padova isochrones give comparable optical/NIR colours 
V-I and V-K, but are offset from the updated Padova94 
models by ~ -O.lmag (V-I) and ~ -0.6mag (V-K), in the 
sense that the updated Padova94 models are redder. 

— The BASTI "non-canonical models" (i.e. with core con- 
vective overshooting during the H-burning phase) are 
closer to the updated Padova94 isochrones than their 
"canonical models" (i.e. without overshooting). 

— The mass lost due to stellar evolution differs by 2% (new 
Padova) - 7% (BASTI) when compared to the updated 
Padova94 isochrones. 

— The relative effects induced by the preferential mass loss 
(i.e. the difference between models with and without 
the effects of cluster dissolution) are qualitatively ro- 
bust against the choice of isochrones. Small quantitative 
differences are present. However, they tend to be even 
stronger for the new isochrones than for the updated 
Padova94 isochrones. 

For two example colours (B-V and V-K), the time evolu- 
tion for a standard SSP (i.e. without cluster dissolution) at 
solar metallicity is shown in Fig. [121 The blue colour B-V, 
dominated by hot stars mainly on the main sequence, shows 
good agreement between the investigated isochrones. The 
V-K colour for ages younger than «1.5 Gyr is in strong 
disagreement between all three isochrones, up to 2.8mag 
difference around an age of 300 Myr, showing strong differ- 
ences in the treatment of the AGB phase. For older ages, 
the new Padova agrees well with the BASTI isochrones, but 
both are offset from the upgraded Padova94 by Ri0.6mag. 

To account for the uncertainty in the choice of 
isochrones, we will build grids of dissolving cluster models, 
based on the BASTI, the new Padova and the Berte lli et al.l 
(|2008t ) (once the extension to higher masses is published) 
isochrones, respectively, and release it on our w e bpage . 
These models will also employ the iKalirai et al.l 
initial-final mass relation. 



4.4. Comparison with earlier work 

The new models pre sent ed here supersede our earlier work 
(|Lamers et al.l 120061 ) . In iLamers et al.l (|2006( l we approx- 
imated the changes in the mass function by a time- 
dependent lower mass limit (i.e. assuming that only the 
lowest-mass stars are removed from the cluster, while 
higher-mass stars might only be removed by stellar evo- 
lution) and scaled our models to match the total mass in 
stars with M < 2 Mq with the BM03 simulations. Recently, 
this a pproach has been improved by Kruijsscn & Lamers 
(j2008f ) by incorporating the effects of stellar remnants for 
clusters of different initial masses and different total dis- 
ruption times for a range of metallicities. In that paper, 
the consequences of various physical effects on the photom- 
etry and M/L ratios have been investigated, e.g. initial mass 
segregation, the role of white dwarfs and neutron stars, the 
role of metallicity etc.. 
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Fig. 12. Integrated B-V (left panel) and V-K colour (right 
panel) for a Kroupa IMF, solar metallicity, no cluster 
dissolution and 3 different sets of isochrones: the "up- 
graded Padova94" (our standard models, black lines), the 
"BASTI" isochrones (red fines, iPietrinferni et al.l WM) 
and the "new P adova" isochrones (blue dotted lines, 
iMarigo et a"Lll2008f ). 



Due to the normalisation procedure, the total masses of 
the earlier models differs negligibly from the new models. 
However, the number of bright stars in th e new models 



decrea ses slower than in the old models by Lamers et all 
(|2006f l (the models bv lKruijssen fc Lamersll2008l represent 
already an improvement over the earlier work, and are more 
consistent with the work presented here). Hence the new 
models are brighter than the old models, especially for short 
disruption times. Consequently, the new mass-to-light ratio 
is lower, by 20 - 40%. 

The new models are redder than the older ones, because 
the changing slope of the mass function slowly depopulates 
the (blue) main-sequence turn-off region already early on. 
Contrary, in the older models stars in the main-sequence 
turn-off region are removed more abruptly when the lower 
mass limit reaches the turn-off mass. 

Before the lower mass limit reaches the turn-off mass in 
the old models, colors get slightly bluer for a short time, be- 
cause almost all main sequence stars redder than the turn- 
off have been removed by then. This feature is not present 
in the new models, due to the more gradual mass loss. In 
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addition, the old models show a strong reddening in their 
final stages, as the star cluster contains exclusively red gi- 
ants/ AGB stars (plus stellar remnants). Also this feature is 
not strongly present in the new models, as the mass func- 
tion even close to total disruption covers a wider range. 



5. Conclusions 

We presented a novel suite of evolutionary synthesis mod- 
els, which accounts for the dynamical evolution of star clus- 
ters in a tidal field in a realistic manner. The dynami- 
cally induced changes in the stellar MF within the cluster 
and the overall mass loss of stars from the cluster into the 
surrounding field population is consistently taken into ac- 
count. The models are made publicly available on our web- 
pages http:/ /www. phys.uu.nl/~anders/data/SSP_varMF/ 
and http://data.galev.org for general use. 

Based on the simulations presented in BM03, we im- 
proved the parametrisation of the time evolution of the MF 
slope changes. We then combined this new description of 
the MF slopes with our GALEV evolutionary synthesis mod- 
els. The resulting models, calculated for a range in mctal- 
licities and total cluster disruption times, were shown to 
deviate significantly from the canonical evolutionary syn- 
thesis models which neglect the effects of dynamical cluster 
evolution. Depending on the total cluster disruption time 
and the colour index under investigation, differences up to 
0.7mag (and in a large number of cases exceeding O.lmag) 
were found. These deviation were shown to lead to signifi- 
cant misinterpretations of observation. E.g. cluster age de- 
terminations can be wrong by 20 - 50%, in extreme cases 
by up to a factor ~2 - 3. These deviations were found to 
depend strongly on the filter combination used to derive the 
ages: combinations including near-IR filters tend to be more 
sensitive to the changing MF, while for large wavelength 
coverage and/or large numbers of filters the deviations are 
still significant but generally smaller. 

Also the M/L ratios are strongly affected, and therefore 
photometric cluster masses derived from observations. For 
the largest part of a cluster's lifetime the M/L ratios are 
significantly below the canonical values (by up to a factor 
~3 - 7). In late stages of cluster dissolution the M/L ratios 
exceed the standard values, as the cluster mass gets in- 
creasingly dominated by stellar remnants. This period can 
last for up to ~16% of the cluster's total disruption time. 
In both cases, the M/L ratios are strongly time-dependent. 
For fixed cluster age and/or fixed local disruption time, the 
dependence of M/L ratios on the presently observed clus- 
ter mass was investigated. They are broadly consistent with 
observations, although the observations show large scatter 
and uncertainties. 

Our results confirm the trends in the evolution of colour 
and mass-to-light ra t ios of diss olving clus t ers, o btained by 
iKruiissen &: LamersI (|2008D and lKruiissenI (|2008D . who used 
a simplified description of the changes in the mass function 
due to the preferential loss of low-mass stars in star clusters. 

While the absolute values of our results depend on our 
choice of input physics, the general behaviour is robust 
against these choices. We will update our models whenever 
better input physics becomes available. 
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